Spin Filtering and Entanglement Swapping through Coherent Evolution of a Single Quantum Dot 
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We exploit the non-dissipative dynamics of a pair of electrons in a large square quantum dot to perform 
singlet-triplet spin measurement through a single charge detection and show how this may be used for entangle- 
ment swapping and teleportation. The method is also used to generate the AKLT ground state, a further resource 
for quantum computation. We justify, and derive analytic results for, an effective charge-spin Hamiltonian which 
is valid over a wide range of parameters and agrees well with exact numerical results of a realistic effective-mass 
model. Our analysis also indicates that the method is robust to choice of dot-size and initialization errors, as 
well as decoherence introduced by the hyperfine interaction. 



Introduction - Realizing quantum information and compu- 
tation tasks in solid state systems, particularly quantum dots 
(QDs), has attracted a lot of interest in recent years. Elec- 
tron spins in QDs are promising candidates for the physical 
implementation of a qubit [1] due to their long coherence 
times J2]. Initialization, manipulation, and readout of elec- 
tron spins have already been demonstrated 

Hi 

and ideas 

exist for quantum gates based on single qubits encoded in two 
QDs 15[]. As it is timely for "proof of principle" demonstra- 
tions of multi-qubit processes, it would be highly desirable to 
establish a coherent two qubit process in a single quantum dot. 

Bell measurement is a key ingredient that makes possible 
some important tasks such as teleportation |0] and entangle- 
ment swapping (7]. In this Letter, we propose a mechanism 
for singlet-triplet measurement based on the coherent dynam- 
ics of two electrons in a large square QD, followed by a single 
charge detection. Such spin-filtering will give a perfect Bell 
measurement in the S z = subspace of two spins. This pro- 
jection is made possible due to the existence of a ground man- 
ifold of two singlets and two triplets, separated from higher- 
lying states by a large energy gap. To a very good approxi- 
mation this enables the low-energy coherent dynamics to be 
confined to the ground manifold in which the singlets rotate 
around the quantum dot whereas the triplets are frozen at their 
initial locations. By initializing the system in an unentangled 
superposition state we are then able to project onto a singlet or 
triplet state simply by a charge measurement to detect whether 
or not the charge has moved during the evolution. We use this 
property to propose some quantum information applications 
such as entanglement swapping and generating the Affleck- 
Kennedy-Lieb-Tasaki (AKLT) state, which is a resource for 
measurement-based quantum computation (H]. 

Recently a dissipative method for singlet-triplet measure- 
ment has been implemented experimentally J3]. In this 
method a double QD is prepared with one electron in each 
QD, and after lowering the barrier one of the electrons will 
hop to the other QD provided that they are in a singlet state. 
As the singlet state is produced by a dissipative decay, there 
is no set time at which the electron will hop and the timescale 
for dissipative relaxation is usually longer than coherent evo- 
lution in the same range of energy. In our coherent mecha- 



nism, however, the operation time is precisely known and the 
filtering takes place well within the spin-coherence time. 

From a practical perspective a large square QD is easier to 
fabricate than a small one and will also be modeled more ac- 
curately by our effective Hamiltonian, since the energy gap 
between the ground manifold and the lowest excited states in- 
creases rapidly with dot size, making the ground manifold in- 
creasingly isolated. On the other hand, as the absolute sizes 
of the singlet-triplet splitting in the ground manifold fall expo- 
nentially with dot size, large QDs have slower operation times 
and are more susceptible to errors. There is thus a trade-off 
between these factors, favoring QDs of intermediate size. Our 
simulations show that for square QDs of L = 200 — 800 nm 
our effective Hamiltonian is sufficiently accurate, and oper- 
ates at frequencies within the range achieved in recent current 
experiments 

Effective Hamiltonian - We consider a system of two elec- 
trons held in a square semiconductor QD with a hard-wall 
boundary, which can be realized in experiment by gating a 
two-dimensional electron gas (2DEG) at a heterojunction in- 
terface. The spectrum of this system is determined by the 
competition between the kinetic energy (~ 1/ L 2 ) of the elec- 
trons and the Coulomb repulsion (~ 1/L) between them. In 
small QDs the kinetic term dominates, and the charge density 
peaks at the center of the dot (like non-interacting particles). 
Conversely in large dots, when the Coulomb interaction dom- 
inates, the energy of the system is minimized by the electrons 
localizing in space to minimize the electrostatic interaction 
energy. In analogy to the concept of the Wigner crystal state in 
bulk two-dimensional systems, these highly-correlated quasi- 
crystalline states are termed "Wigner molecules". 

Assuming an effective mass m* for the electrons the square 
QD is modeled by: 
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[V? + V 2 ] +V(r 1 ) + V(r 2 ) 



47re|ri - r 2 | 
(1) 

where V(r) is the confining potential. We choose this to be 
hard-wall with exact square symmetry, though our results will 
not qualitatively change under small deviations from perfect 
symmetry. The last term in Eq. ([T]i represents the Coulomb re- 
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Figure 1: (Color online) Eigensystem of a GaAs dot with a side- 
length of L — 800 nm, obtained by exact diagonalization of the 
effective-mass Hamiltonian (Eq. Q}. (a) The lowest two multiplets 
of states; singlets are shown with solid (blue) lines, triplets with 
dashed (red) lines. We consider only the dynamics of the lowest mul- 
tiplet, consisting of two singlets (\Si) and 15*2}) with two degenerate 
triplets lying between them, (b) Charge distribution of the ground- 
state, showing the formation of a Wigner molecule, with peaks la- 
beled abed near the dot corners, (c) Charge distribution of the sym- 
metrized singlet state |1) = (|Si) + |52})/\/2> localized about bd. 
(d) Charge distribution of the antisymmetrised singlet state |2), lo- 
calized about ac. 

pulsion between the two electrons. In the strongly-correlated 
regime, in which the size of the square is large compared with 
the effective Bohr radius (~ lOnm in GaAs), eigenstates 
of this simple Hamiltonian are extremely demanding to ob- 
tain exactly. We show in Fig. Q2a) the low-lying energy spec- 
trum of a GaAs QD with side-length 800 nm, obtained by 
diagonalizing the full two-electron Schrtidinger equation. We 
see that two degenerate triplets (|n), n = 3,4, ...,8) sit ap- 
proximately (but not precisely) midway between two singlets 
(|5 1 ( 2 ))), while all these 8 states are separated from the next 
multiplet of eigenstates by a relatively large gap. The charge 
distribution for the ground-state \Si) is shown in Fig. Q2b), 
and clearly shows how the charge density strongly peaks near 
the corners of the QD. One can better appreciate the form of 
the states by defining linear combinations of the two singlets 

11} = (|Si> + |Sa»/V2 = |Sf)|,n (2) 

|2) = (\S 1 )-\S a ))/V2 = \$i)\4,-), 0) 

where = (| fj,) - | |t))/\/2 is the singlet spinor, 

and |$f( 2 )) is the symmetric spatial component of the two- 
electron wave function. In Fig. [He) and [TJd) we plot the 
charge distribution of these states, clearly showing how they 
are localized at diagonally-opposite corners of the QD. For the 
triplets we adopt a similar labeling scheme 

|3) = |^)|^+),|4} = |$ 2 4 }|^ + } ) |5> = |^}|tt} 1 (4) 
|6) = |*£>|tt>,|7> = U)A&) = \*t)\ U), (5) 




Figure 2: (Color online) Gate structure for a large QD (central shaded 
square), connected to two smaller QDs (pink circles) at opposite cor- 
ners. 

where \i/>+) = (| U) + I lt))/v^, and (|$£)) is the 
anti-symmetric charge distribution, which resembles that of 
the states |1) and |2), being peaked at the same sites ac (bd). 
Note that while the triplets \n) (n = 3, 4, 8) are eigenvec- 
tors of H, the singlets |1) and |2) are not. 

We can immediately write down an effective Hamiltonian 
for the low-lying energy eigenstates 

8 

ff eff = -Ai|5i)(5i| + A 2 \S 2 )(S 2 \ + E £ \n)(n\, (6) 

n=l 

where Eq is the energy of the two degenerate triplets, and Ai 
(A2) is the energy separation between the triplets and \Si) 
(I1S2)). By restricting ourselves to the ground manifold, and 
using the sum rule Yl n =i \ n )( n \ = I> tne effective Hamilto- 
nian may be written in the charge-spin form 

H eS = E I - A(|l)(2| + |2)(1|) + J( Sl • s 2 - 1/4) , (7) 

where J = (A 2 - Ai)/2 and A = (Ai + A 2 )/2. This 
has the simple physical interpretation that the Coulomb re- 
pulsion pushes the electrons to diagonally-opposite corners, 
giving two charge states for each combination of spin. Whilst 
in the corners the spins of these electrons have an effective 
Heisenberg exchange interaction, with exchange constant J, 
and they may tunnel from one charge state to the other with 
amplitude A. 

Dynamics - We now consider the time evolution of two 
electrons injected into the square dot so that one is located 
near a and the other near c (as labeled in Fig. [lib)). This 
could be achieved in principle using surface gates as shown 
schematically in Fig. [2] Initially an electron is localized in 
each of the small dots adjacent to the large dots. These elec- 
trons are then transferred to the large dot by lowering barriers 
using gates Gi,Gs and Gi,G§ which are subsequently re- 
stored to their previous potentials after electron transfer has 
completed. If both electrons have the same spin, i.e. total 
S z = ±1, then this spin will not subsequently change with 
time under the coherent evolution of the Hamiltonian ^} and 
the two electrons will therefore remain close to their parent 



corners. However, if the two injected electrons are of oppo- 
site spin, then the state after injection will be an equal super- 
position of a singlet state and an S z = triplet state, which 
will subsequently change with time. Let us consider the state 
in which a spin-up electron is injected at corner a and a spin- 
down electron at corner c. We may approximate this state by 
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Note that both components correspond to spin-up at a and spin 
down at c since $f + $^ ~ except when ri ~ r a , r 2 ~ r c 
and $f — <3?^ ~ except when ri ~ r c , r 2 ~ r a . Hence this 
state is unentangled. Under the Hamiltonian ©, the time- 
evolution of \ip(0)) can be determined analytically as 
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V2 



e lJt (cos(Ai)|l) +isin(Ai)|2)) 



|3)], 
(9) 

1, We see directly that at time t* = 
1, we have a superposition of 



choosing units with h 
tt/2 A, for which sin(i*A) 
the two states \2) and |3) with the same probability of finding 
either of them. At time t*, therefore, a simple single charge 
detection at any corner (let us say b) will project \ip(t*)) into a 
singlet (with the electrons in corners b and d) or a triplet (with 
electrons remaining in corners a and c). Hence, if we project 
the state into a singlet then it oscillates between corners bd and 
ac. Conversely, if we project it to the triplet then it is frozen 
in the corners a and c. 

The probability of detecting the singlet state at time i, start- 
ing in the S z = subspace, is P 2 = \(2\i'(t))\ 2 = \ sin 2 At. 
Thus P 2 oscillates harmonically with maximum probability 
1/2 but independently of the exchange, J, which simply in- 
duces a phase factor in the singlet component of the wave 
function. This independence of J implies that our method 
of 'filtering' the singlet by measurement is robust to the size 
of the dot, for which the ratio J/ A falls exponentially with in- 
creasing dot size lfToh . This is not the case for other overlaps. 
For example, the probability of finding the initial state is 



P m = \(i>(o)\m)\ 2 = 



1 + cos 2 At + 2 cos Jt cos At 



(10) 
0) does the 



which shows that only for special cases (e.g. J 
system return to its starting state. 

Applications - The ability to make singlet-triplet measure- 
ments paves the way to implement some quantum computa- 
tion tasks such as entanglement swapping, or equivalently, 
teleportation. To achieve these we generate two singlet pairs 
outside a square QD as shown in Fig. |3 a )- These pairs may 
be generated via surface gates in a similar fashion to those 
shown in Fig. [2] in which electrons are transferred from the 
surrounding 2DEG reservoir. The singlets are formed sim- 
ply by cooling the system U We then push one electron 
from each singlet pair into the big square QD as shown in Fig. 




(b) 




(d) 
►•1 




Figure 3: (Color online) (a) two small QDs, with a singlet pair in 
each, beside a large square QD (dashed lines denote entanglement); 
(b) One electron from each singlet is pushed into the square QD; (c) 
Entanglement swapping; (d) Scaling up the system to an array of 
QDs. 



[3b). We now have two electrons in the corners a and c in 
the square QD and after time t* we measure the charge at one 
corner. With probability of 1 /4, the state of the electrons in 
the square QD collapses to a singlet at sites bd. In this case 
the two remaining electrons in the small QDs get entangled as 
another singlet, as shown in Fig. f3 c )- This process is called 
entanglement swapping (or the teleportation of entanglement) 
and generates entanglement between distant particles. This 
scheme can be scaled up through a geometry shown in Fig. 
0d) where a series of empty square QDs are arranged between 
small QDs containing electron singlet pairs. By pushing one 
electron from each small QD to its neighboring square QD, 
one makes all small QDs empty except the two which termi- 
nate the array, that each hold one electron. Dynamical singlet- 
triplet measurement on all the square QDs generates a singlet 
between the electrons held in the terminating small dots when 
the result of all measurements is singlet. The probability of 
having this is (1/4) w , where N is the number of square QDs. 

When the result of measurement in Fig. [3b) is a triplet, 
rather than a singlet, we can generate the so-called AKLT state 
1 1 111 . Originally this was introduced as the ground state of 
the AKLT Hamiltonian ll ill , which models the interaction of 
a series of spin-1 particles with two spin- 1/2 particles at the 
boundaries of a chain. The AKLT ground state can be gen- 
erated by again starting with a series of spin- 1/2 singlets in 
small QDs but this time, projecting two particles of neigh- 
boring singlets into a triplet to represent their spin-1 nature. 
This occurs with probability 3/4 when the result of the mea- 
surement in Fig. [3b) is a triplet. This can also be scaled 
up with the geometry shown in Fig. [3d), with probability of 
success is (3/4) w that all square QD states will be in a triplet 
state. The AKLT state can be used as resource for ground- 
code measurement-based quantum computation (H]. 

Gate Errors - The above results for the time-development 
of the initial state are exact, requiring only the energy pa- 
rameters J and A, which can be obtained directly from the 
eigenenergies of the ground manifold of the effective-mass 
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L (nm) 


A (meV) 


J (meV) 


|a| 2 


n 2 


EhffaeV) 


100 


0.814 


-0.243 


0.441 


5.23xl0" 2 


1.74 


200 


0.145 


-4.363xl0~ 2 


0.445 


3.63xl0~ 5 


7.76 xlO -1 


400 


2.11xl0~ 2 


-5.05xl0~ 3 


0.420 


1.21 xl0~ 3 


3.88X10- 1 


800 


2.08xl0 -3 


-2.20xl0" 4 


0.453 


2.78xl0 -4 


1.94xl0 _1 


1600 


9.34xl0" 5 


-1.66xl0~ 6 


0.490 


6.02xl0~ 6 


9.69xl0~ 2 



Table I: Physical parameters for a GaAs QD. \a\ 2 and |/3| 2 (Eq. (TT\ ) 
are the projection of the initial state onto the singlet states |1) and |2) 
by applying a gating potential of 0. 1 V. 



Hamiltonian Eq. (HJ. However, these results are somewhat 
contrived in that the starting state lies precisely within the 
Hilbert space of the ground-manifold, and therefore remains 
within this ground manifold under time evolution. In any re- 
alistic situation these conditions will not be met and in par- 
ticular the starting state will deviate from the idealized form, 
Eq. (8). It will contain small admixtures of the other base 
states in the ground manifold and excited singlet states. These 
admixtures will increase with decreasing dot size but should 
still give small errors for L > IOob, say. We can derive ex- 
pressions for the fidelity starting with a more realistic state, 
|V>(0)). This could be produced, for example, by applying a 
positive potential to gates located near the sites a and c. In 
the numerical calculations, this was modeled by dividing the 
square dot into four quadrants and applying a constant pos- 
itive potential to the two diagonally opposite quadrants that 
contain the corners a and c. In this scheme setting the gating 
potential to 0.1 V yields values for the overlap (ip(0)\ip(0)) 
of 0.80, 0.940, and 0.97 for QDs of L = 200 nm, 800 nm 
and 1200 nm respectively, which are reasonably close to unity, 
and could be enhanced further by using more elaborate gating 
potentials. We may derive an expression for the fidelity with 
this more realistic initial state by expanding | V'(O)) in terms of 
\ip(0)), (|1) — |3})v2 and the remaining eigenstates of the full 
effective-mass Hamiltonian. After time evolution and projec- 
tion onto 12) we obtain 



P| = |(2|^(t))| 2 = (a sin At) 5 



2a/3 sin Jt sin At + f3 2 

where a = (1^(0)) and f3 = (2|^(0)). Note that Pf is 
independent of excited states, and since \a\ 2 ~ 1/2, |/3| 2 ~ 0, 
it is robust to gate errors. This is illustrated in Table I where 
we see only small deviations from the ideal P2, even for the 
smallest dot of L = lOOnm, the main effect being a suppres- 
sion of the maxima and enhancement of the minima. 

Charge measurement - For simplicity we have so far as- 
sumed that charge detection may be made on a timescale much 
less than the coherent charge evolution time t*. Typical val- 
ues of t*, however, being of the order of nanoseconds for our 
parameters (see Table I) are challenging to measure directly 
in experiment. For practical implementation, we propose a 
similar scheme to Ref. H, which is able to achieve an ac- 
ceptable time resolution. At the moment of measurement we 
restore the quadrant gate-potentials (used previously to ini- 



tialize the system) to freeze the dynamics of the electrons. A 
strong charge measurement at one of the corners of the QD 
can then be made to project the state into a singlet or triplet. 

Charge dephasing - Charge dephasing reduces the coher- 
ence between |1) and |2) in Eq. (0, but since our measure- 
ment projects onto these states anyway, it does not funda- 
mentally affect our scheme. By damping the sinusoidal os- 
cillations between |1) and |2), charge dephasing also reduces 
P 2 {t*) = \{2\4>{t*)}\ 2 such that in the extreme case of very 
strong decoherence it goes to 1/4. In this case if |2) is detected 
successfully the scheme is completed as before, giving entan- 
glement swapping. Otherwise, we end up with a superposition 
of |1) and |3), as in the initial state (except for the amplitude 
of 1 1) being reduced), which again undergoes damped oscilla- 
tions. By repeating this process one can reliably (with expo- 
nential improvement according to number of trials) discrimi- 
nate between singlets and triplets in the initial state. However, 
due to our fast dynamics this extreme case is very unlikely. 
As an example, for L = 400 nm we have t* = 0.3 ns, which 
is safely below the dephasing time T2 ~ 1 — 2 ns in a system 
with comparable size Jgt]. 

Hyperfine interaction - Another source of decoherence is 



the hyperfine interaction between electrons and nuclei 11211 . 
This can be estimated by replacing the effect of the nuclei with 
an effective magnetic field B* coupled to the electron spin as 
where j e = gps/fi and a = (a x ,a y ,a z ) 
are the Pauli matrices. B* has a Gaussian random distribution 
with a variance B n iflill . The major effect of the hyperfine 
interaction is to mix the spin-subspaces. Due to the fast evo- 
lution, the first maximum of p2(i) is relatively unaffected for 
typical energy values of Ehf = hj e B n , given in Table I. As 
an example, for L = 400 nm (t* = 0.3ns), the decoherence 
time scale of hyperfine interaction is 10.7ns. 

Conclusions - We have shown that the dynamics of a pair 
of electrons in a large square QD can be used to perform 
singlet-triplet spin measurement using just a single charge de- 
tection. This is accessible to current technology and unlike 
previous schemes, it is fast, deterministic and coherent. Re- 
peating the singlet-triplet measurement to a chain enables en- 
tanglement swap over a distance and the generation of AKLT 
state in a way that would enable proof of principle quantum 
information experiments. Furthermore, coherent evolution of 
the system is considerably faster than the dephasing time T2 
imposed by the hyperfine interaction. Our low-energy analytic 
description is valid for a wide range of parameters, particu- 
larly for typical experimental values of the QD parameters. 
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